Fair sharing of resources in a supply network with constraints 
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This paper investigates the effect of network topology on the fair allocation of network resources 
among a set of agents, an all-important issue for the efficiency of transportation networks all around 
us. We analyse a generic mechanism that distributes network capacity fairly among existing flow 
demands. The problem can be solved by semi-analytical methods on a nearest neighbour graph with 
one source and sink pair, when transport occurs over shortest paths. For this setup, we uncover 
a broad range of patterns of intersecting shortest paths as a function of the distance between the 
source and the sink. When the number of intersections is the maximum and the distance between 
the source and the sink is large, we find that a fair allocation implies a decrease of at least 50% from 
the maximum throughput. We also find that the histogram of the flow allocations assigned to the 
agents decays as a power-law with exponent -1. Our semi-analytical framework suggests possible 
explanations for the well-known reduction of the throughput in fair allocations. It also suggests 
that the combination of network topology and routing rules can lead to highly uneven (but fair) 
distributions of resources, a remark of caution to network designers. 
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I. INTRODUCTION 

Transportation networks carry the flow of sonic com- 
modity or entity using routing rules that aim to be effi- 
cient and avoid congestion [J Recently, researchers 
have focused on the design of optimal transportation net- 
works , the characterisation of robustness and dam- 
age in transport networks LDtSI, and routing strategies 
that minimize congestion [10l4l4l |. Despite these efforts, 
little is known about the effect of network topology and 
various routing protocols on the competing interests of 
network users. Indeed, transport frequently involves a 
very large number of agents, and the operation of the 
network has to be managed in a smooth and fair way. 
Moreover, a socially efficient solution, i.e., one that max- 
imizes the sum of the utilities of individual agents, can 
be difficult to implement because it may be perceived as 
unfair to some of the agents 

The fair division of network capacity among agents is 
a problem related to the procedure of dividing a cake 
fairly. Cake cutting, or the fair division of a divisible 
good among a set of agents, can be analysed u sing several 
mathematical interpretations of 'fairness' |19l - [23 |. When 
all N agents agree on the value of each portion, the prob- 
lem becomes trivial and each agent receives a piece of size 
1/iV. If the transport routes intersect along one edge 
only, the problem of dividing the edge capacity is akin 
to cutting a cake. However, in most situations of in- 
terest, distinct agents share transport routes in sophisti- 
cated ways, and the problem is then how to allocate fairly 



the available capacity on the intricate web of congested 
edges. Two prominent strategies to address this conun- 
drum are the Max-Min Fairness (23j and the Proportional 
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Fairness |2J| allocations, and this paper focuses on the 
former. A set of flows is Max-Min Fair if no flow can 
be increased without simultaneously decreasing another 
flow that is already less than or equal to the former. To 
oversimplify, an allocation is Max-Min Fair if the wealthy 
can only get wealthier by making the poor even poorer. 

The related problem of network congestion control 
aims at devising mechanisms for the allocation of re- 
sources in transport systems [l^, HH . There is a long line 
of work on addressing the challenge of congestion control 
with the Max-Min Fair allocation, mainly by researchers 
working on communication networks. Indeed, the fair al- 
location of capacity has been extensively studied in the 
context of flow control in IP packet switchin g n etworks, 
such as the Internet and ATM networks [23l - |28| . More- 
over, researchers have taken a wide range of approaches 
to the allocation of communication rates, encompassing 
both non pre-allocated [29l - l3]| and fixed [2^, [s^l paths. 
Despite these efforts, to the best of our knowledge, pre- 
vious approaches have focused on the design of algo- 
rithms In [13, [i^-dS [H-il, or the computation of 
upper bounds of network throughput [l8| . In contrast, 
our results offer a proof of concept that flow control can 
be analysed semi-analytically in large networks, at least 
regular ones, and reveal insights into the interplay of net- 
work topology, capacity, route selection and distance be- 
tween a source and a sink, and how these elements in 
turn affect congestion in fair networks. 

The choices of network topology and routing place ob- 
vious constraints on the fair allocation of network re- 
sources. If the routing is such that paths do not intersect. 
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each path is assigned the minimum capacity among all 
its edges and the throughput is the max- flow. In general, 
however, the choice of routing has an impact on the num- 
ber of path intersections and thus on the reduction of net- 
work throughput in fair allocations. One possible choice 
of routing is to allow transport over all possible s-t paths. 
In such case, the Max-Min Fairness algorithm is typically 
solved by an iterative LP model ^30,, , and the fair allo- 
cation to one source-sink pair gives the same throughput 
as max-flowfs^- However, the model is cornputationally 
heavy when the number of s-t pairs is large [36[ . Another 
possibility, followed in the vast majority of work on Max- 
Min Fairness [H, [H [Ulli, 113, Isl] and in this paper, 
is to fix the paths off-line before the algorithm is run. 
Fixing the paths is justified, although network managers 
arc confronted with dynamic real-time traffic conditions, 
because it is frequently impractical to update the rout- 
ing schemes according to changes in the distribution of 
traffic [1^ Hence, network operators often design 

the best possible static network conditions for a set of 
estimated average or worst-case demands [11] ■ Here we 
choose to analyse routing over shortest paths, because 
this is a highly schematic way to model several transport 
processes in the real world, which is frequently used in 
the literature [H, H^] . Some real- world examples moti- 
vating its usage are packet transfer on the Internet 4l|. 
spatial distribution networks such as sewer systems 4^ , 



gas pipeline networks 
networks [4311. 



commuter rail j42l | and subway 



II. FAIR DIVISION ON NETWORKS 



A. How to share network capacity fairly in supply 
networks 




FIG. 1. (Colour online). Illustration of the Max-Min Fair 
(MMF) allocation on a small water pipeline network with 
three source nodes (identified by water taps), two sink nodes 
(identified by buckets), and pipelines with unit capacity, (a) 
Shows the initial configuration, before the application of the 
MMF algorithm. Five paths, connect the sources (taps) to 
the sinks (buckets), (b) At iteration one, the flow is in- 
creased uniformly on all taps, until the capacity is reached 
on pipeline 61,4. The path flows that cross ei,4 are frozen, 
the corresponding taps can no longer add to network flow 
and are thus removed from the flgure, and we say that the 
pipeline is saturated (saturated pipelines are marked in dark 
gray colour), (c) and (d) The process is repeated for the un- 
saturated pipelines, until pipelines 64,5 and 63,5 are saturated 
at iterations two and three, respectively. The sink inflow /4 
and /s is displayed on each panel. The flnal allocation on 
panel (d) is Max-Min Fair, because no path flow can then 
be increased without decreasing another path flow already 
smaller. 



The Max-Min Fair (MMF) allocation can be found 



with the water-filling algorithm |37| . The algorithm 
starts with all path fiows from zero. Path fiows are then 
increased equally for all paths, until one or more edges 
arc saturated. We say that such edges are bottlenecks, 
and paths that pass through bottlenecks are also said to 
be saturated. The path flow of saturated paths is the 
Max-Min Fair flow allocation of the path. The algorithm 
is repeated until all paths are saturated, that is until each 
path passes through at least one bottleneck edge. 

We give a visual representation of the water-filling al- 
gorithm in Fig. [ll where a water network connects taps 
(sources) to buckets (sinks) with individual paths on 
pipelines of unit capacity. Here, the algorithm proceeds 
by simultaneously opening the taps at the sources (see 
Fig. [1^). In the first iteration of the water-filling algo- 
rithm all paths get an equal flow of 1/3 and the edge 61^4 
becomes saturated (see Fig. [T}j) . In the second iteration 
of the algorithm, the tap at node 1 cannot be opened 
any further, as the unit capacity of 61^4 has been fully 
used. Nevertheless, the taps at nodes 2 and 3 can still 
be opened further, and the path flows of paths passing 



through 62,4 and 63,5 are increased to 2/3 (see Fig. [T};). 
After the second iteration, the tap at node 2 cannot be 
opened further as the edge 64,5 is saturated. On the third 
iteration, the tap at node 3 can be opened further, until 
the path flow on the path passing through 63,5 is 1 (see 
Fig. [TJl). The algorithm then stops and the allocation 
is Max-Min Fair because no path flow can be increased 
without decreasing another path flow less or equal to it. 
The analogy with cake-cutting can now be established: 
at each iteration of the algorithm, the available capacity 
on the bottleneck edge is shared equally among the paths 
that cross it and that have not been allocated a share of 
any edge capacity in previous iterations. In other words, 
allocating the network capacity can be seen as cutting 
a sequence of 'cakes' (the capacity of bottleneck edges), 
and the Max-Min Fair allocation establishes the sequence 
of 'cakes' to be shared equally among the paths that cross 
the bottleneck edges. We are now ready to formalise the 
iterative MMF algorithm [lllli,!!^. To do this, we first 
need some mathematical definitions. 
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B. The Max-Min Fair algorithm explained formally 

Let G — (y, c, 5, T) be a connected and undirected 
network with node-set V, edge-set E, an edge capacity 
function c : E ^ R J , and a set of source and sink pairs 
{si,ti) G S xT with i = 1, • • ■ ,M. Each source and 
sink pair {si,ti) is connected by a set of Si ~ U paths 
Pi = {P{i,i),--- ,P(i,fci)}, where P = lifi^Pi is the set 
of all source to sink paths on the network. The undi- 
rected edge Cij £ E connects nodes i and j, and we will 
omit the subscript indexes when we are not referring to 
a specific edge. All edges e G E{p{^i^k)) of a s; — U path 
P(i,k) (1 ^ * ^ M and 1 < fc < fci) transport the same 
path flow fpi^. G M.Q between s,; and ti. This implies 
conservation of the flow at the nodes. Different paths 
can share an edge, even to perform transport in differ- 
ent directions {e.g., during distinct time intervals) psj . 
The set of paths passing through edge e & E is P{e) = 
{P(z,k) e -P : e e £:(p(,_fc))}. The flow on edge e £ E is 
the sum of the path flows on paths passing through the 
edge, and is given by F{e) = Ep(.,,)eP («^)/p(..»=) ' 
where Sp^. ^^^{e) is 1 if e S E{p(^i^i.-^) and is otherwise. 
A flow is then feasible if < F{e) < c(e) for all e & E, 
where c(e) is the capacity of edge e. Formally, a set {/} 
of path flows on P is Max-Min Fair, if it is feasible and if 
for any other feasible set {/'} of path flows on P, there 
exists a path G P : fp^^ ^^ > fp^^ ^^ implies that 

there exists another path P(jj) E P '. fp^ < fp^ - and 

fpu.i) — fpi^.k) ■ 

We define P*^*-* to be the set of paths on the network 
at iteration i of the MMF algorithm, and P'-'^(e) to be 
the subset of paths in P^^'> that pass through edge e. 
Before we start the algorithm, we assign P(i) ^ P and 



.(1), 



e) = c(e) for all e G E, and a path flow /, 



(0) 



to 



each path P(j^k) & P*-^-*- We then start the algorithm and 
initialize the iteration counter i — I. 

In the first step of the MMF algorithm, for each edge 
e with non-zero capacity that belongs to at least one 
path, we define the edge capacity divided equally among 
all paths that cross the edge at iteration i of the MMF 
algorithm as 



,/,W(e) = c«(e)/ P«(e) 



(1) 



for all c(')(e) =^ 0. We then find the minimum of (/)'^*'(e), 
given by 

A/(0 



min 



(2) 



In the second step of the MMF algorithm, we increase all 
path flows of paths in P^*^ by A/*^*', such that 



f(i) 
■'P(j.k) 



i^;7^)-fA/W ifp(,,.)GPW 



M-l) 

JPU.k) 



ifpO-,fc) GP\P(^) 



(3) 



(i) 



The effect is to saturate the set of bottleneck edges E^ 
{eB&E: Ep,,,,,6P(.)^P„,.,(eB)A/W = c«(eB)}, and 



consequently also to saturate the set of paths that contain 
at least one bottleneck edge. Next, we create a residual 
network, by subtracting the capacity used by the path 
flows. 



PU.k)(^P^'^ 



(e)A/ 



ii) 



(4) 



Note that all edges es G Eg will be saturated, that is 
each will have c^*'^^^(eB) = after this step. We also say 
that all paths that contain at least one edge G E^j^^ 
are saturated paths, to mean that their path flow will 
not be increased in subsequent iterations of the MMF 
algorithm. Next, we remove the set of saturated paths 
from P, that is 



pii+l) _ U 



P^'Kes). 



(5) 



We say that p(*+i) is the set of augmenting paths 
because the path flows of paths in p(*+i) can still be 
increased in subsequent iterations of the algorithm. If 
jg Qjj^pty^ increase the iteration counter, 
i ^ i + 1, and we go back to the flrst step, otherwise we 
stop and store the value of i as i*. The Max-Min Fair 
flow on edge e is then the sum of path flows over all paths 
that cross the edge after the algorithm terminates: 



Fie) 



PU.k) ■ 



(6) 



At each iteration, the MMF algorithm increases all path 
flows equally. Thus, the algorithm flnds the bottleneck 
edges, which determine the maximum path flow possible 
for path flows that are the smallest. In other words, 
by finding the minimum of Eq. ([T]) at each iteration, we 
ensure that the smallest path flows are maximized. 



III. MAX-MIN FAIR FLOWS IN NEAREST 
NEIGHBOUR NETWORKS 

A. The network constraints 

Recall from Eq. ([T]) that the capacity available on the 
bottleneck edges at each iteration of the MMF algorithm 
is shared by the paths passing through those edges. Thus, 
on a network with uniform capacity, the MMF flow allo- 
cation depends exclusively on the number of paths pass- 
ing through each edge. To gain a better insight into the 
MMF allocation of flows, we choose to study if-nearest 
neighbour networks where each edge has constant capac- 
ity before the MMF algorithm is applied, that is. 



c(e) 



.(1) 



(e) 



for all e £ E. 



(7) 



The advantage of studying these regular networks with 
constant capacity is that we can flnd simple closed form 
expressions for counting paths. In these regular graphs of 
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FIG. 2. (Colour online). Two layouts of 6- nearest neighbour 
networks. Nodes are numbered by an index from west to east, 
and the distance d that separates two nodes is the difference 
between their indexes. A source s and a sink t pair are placed 
(a) d = 7 and (b) d = 9 nodes apart, respectively. The 
shortest paths between s and t are identified by numeric labels 
on the edges. The shortest path length between the two nodes 
is 1/ = 3 in both panels. There are six s-t shortest paths in (a) , 
but only one in (b). Edges have unit capacity, edge thickness 
is proportional to the MMF flow allocation, and saturated 
edges (or bottlenecks) are drawn in black. Nodes are coloured 
according to their shortest path length from the source. 



even degree K, the nodes are placed on a one dimensional 
lattice with periodic boundary conditions (the lattice is 
a ring). Nodes are then connected to their nearest, next- 
nearest neighbours and so on, up to the constant range 
K/2. A'-ncarcst neighbour graphs were recently the ob- 
ject of intense study as prototypes of networks where 
small amounts of rewiring can generate small- world net- 
works [i^ll^l. To simplify the analysis, we also introduce 
the constraint that transport only occurs along the set of 
shortest paths connecting the source and sink nodes. 



B. Shortest paths in if- nearest neighbour networks 
with one (s,t) pair 

We place a source and sink pair {s,t) on a if-nearest 
neighbour network with \V\ nodes. Nodes are numbered 
by an integer index, which indicates their relative posi- 
tion d < l\V\/2\ from the source. Hence, there are two 
ways of measuring the distance between the source s and 
the sink t. We can either measure the s-t distance as the 
difference d between the indexes of s and t on the ring, 
or as the shortest path distance L between s and t. Since 
each node has K neighbours, there are K/2 nodes at the 
shortest path distance L from the source, each of them 
at a distinct distance d from the soturce. Fig. [2] is an il- 
lustration of the diversity of source to sink (s-t) shortest 
paths when K ~ 6 and L = 3, ion d = 7 and d = 9. 

We can now analyse the smallest distance d between s 
and t for a given shortest path length L. This distance 




1=0 1=1 1=2 ... 1=L-1 1=L 



FIG. 3. (Colour online), (a) Layout of a if- nearest neigh- 
bour network where the s-t pair is separated by ci™'"(i(', L) 
nodes on the ring for a given shortest path length L. Nodes 
are coloured according to their shortest path length from the 
source. Dashed edges connect two nodes at the same shortest 
path length from the source, (b) Rectangular layout of the 
network shown in (a), where all nodes except the source and 
sink are rearranged onto a rectangular lattice K/2 x (L — 1) 
according to their distance from the source. Nodes on this 
rectangular lattice are identified by their {k, I) coordinates, 
(c) A"-nearest neighbour network, where s-t are separated by 
d = d'"'"(is:,L)-hi nodes on the ring (i = 1, . . .,K/2~1), has 
the same structure of shortest paths as a less dense {K — 2i)- 
nearest neighbour network, where the two nodes are at the 
smaller distance d = d™'"(is' — 2i,L). This panel illustrates 
the case i — 1, for which the gray area shows that s-t shortest 
paths no longer cross row k — 1. 

is given by 

d'^'"^iK, d)^[d- {{d - 1) (mod K/2))] (8) 
d'^'"'{K,L) = {L-l)*K/2 + l, (9) 

as is illustrated in Fig. [3|i. To simplify the counting, 
we rearrange the network layout as in Fig. [3]3, where 
each node (except for the source and sink) is identified 
by its row and column coordinates {k,l), 1 < k < K/2 
and 1 < / < L — 1. A s-t shortest path can be vi- 
sualised on such rectangular subgraphs as a sequence of 
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nodes s,W(fei,i),W(fc2,2), ■ • ■«(fci_i,L-i),i, where the nodes 
''^{ki,i) S V arc such that there is an edge between U(fci,j;) 
and vi^ki+i,i+i}, and fci < fc2 < • • • < fc^-i is a non- 
decreasing sequence of rows in the rectangular subgraph. 
In other words, shortest paths contain only horizontal 
(east) and diagonal (south east) edges. The number of s-t 
shortest paths when the pair is at the distance d"""(A', L) 
is the number of ways we can distribute L — 1 identical 
edges to K/2 distinguishable rows. Such allocation can 
be made in 



N{K/2,L) = 



K/2 
L- 1 



K/2 + L-2 
L- 1 



(10) 



different ways, where ((J^)) = ("'*'™~^) is the number of 
ways that m indistinguishable balls can be assigned to 
n distinguishable urns (see the twelvefold way of combi- 
natorics in [H, pp. 31-38]). Equation ((TO)) implies that 
the number of s-t shortest paths between a pair of nodes 
at distance d™^^{K,L) is given by an entry in Pascal's 
triangle as visualised in Fig. |31 



{n(ii' = 2,d)}d=i,...,L|v|/2j = 1 



3 4 5 



{n(A- = 4,d)}d=i,...,L|v|/2j = 1.1, 21,1, .3,1, 4,1, 5.1 ... 

L=\ 2 3 4 5 

{n(A' = 6,d)}d=i,....L|v|/2j = 1,1.1,^,2, 1,6. .3,1, 10 ,4,1, 

L=l 2 3 4 

15,5,1,... 

5 

{n{K = 8, rf)}d=i,...,L|v|/2j = 1, 1, 1. 1,|, 3, 2, 1,10, 6, 3, 1, 
20, 10 ,4,1,35.15,5,1,... 



{n{K = 10, d)}rf=i,...,uv|/2j = 1, 1, 1, 1, 1,|, 4, 3, 2, 1, 

' ' . ' 

L=l 2 

10, 6, 3, 1, 35, 20, 10, 4, 1, 70, 35, 15, 5, 1, . . . 




FIG. 4. (Colour online). The entries in Pascal's triangle are 
the number of s-t shortest paths between a source and sink 
nodes at distance ci™'"(Jk', L) from each other on a A"-nearest 
neighbour network, for a chosen shortest path length L be- 
tween the two nodes. Entries of Pascal's triangle coloured in 
blue illustrate that (g) is the sum of the number of s-t shortest 
paths over the hockey stick. 

We are now able to determine the number of s-t short- 
est paths when the pair is at distance <F^™{K,L). Nev- 
ertheless, this is not the general case. To illustrate 
the point, the number of s-t shortest paths in nearest 
neighbour lattices is represented in Fig. [5] as a sequence 
{n[K, fi)}^;^! |^|y|/2j : which is in itself a batch of mono- 
tonically decreasing sequences. Each integer at a position 
d of the sequence is the number of s-t shortest paths when 
the source and sink nodes are at distance d on the ring. 
Each of the sub-sequences is characterised by a {K, L) 
pair, but so far we have only determined the first element 



FIG. 5. (Colour online). Batch of sequences denoting the 
number of shortest paths for a single {s,t) pair positioned 
on a AT-nearest neighbour network at distance d from each 
other. K is the node degree and L is the length of a shortest 
path between s and t. The position of the number of shortest 
paths for each K value, is the s-t distance d, e.g. there are 
4 s-t shortest paths for K = 8 a.t a. s-t distance d = 5 and 
shortest path distance L — 2 (see node |4 ). The boxes in 
blue and gray show that the structure of s-t shortest paths 
for given K and d can be recovered as a function of less dense 
K-nearest neighbour networks, where the s-t pair is closer (see 
also Fig. [4|. 



of such sub-sequences, which is N{K/2, d™'"). These first 
elements of each sub-sequence follow a hockey stick pat- 
tern in Pascal's triangle, which is illustrated in Figs. [4] 
and [5] for K — 10 and L — 3. Such a pattern tells us that 
the first element of each subsequence can be found from 
the sum of previous first elements for a smaller shortest 
path length L- 1, that is N{K/2, L) = J^fl^ N{i,L~l). 

To understand the structure of the sub-sequences, ob- 
serve that the set of shortest paths on a if-nearest neigh- 
bour network when s and t are placed d = d™™(if, L)-\-i 
nodes apart is equivalent to the set oi s — t shortest paths 
on a less dense (isT — 2i)-nearest neighbour network where 
s and t are placed closer to each other, at the distance 
d™(i^ - 2i,L). For example, if if = 8 and L = 4, and 
the s-t nodes are at distance (i™™(8,4) = 13, increasing 
the distance to c? = 14 docs not change the shortest path 
length L, but the set of shortest paths between the two 
nodes is the same as when the s-t pair are placed on a 
6-nearest neighbour network at distance c?"""(6,4) = 10. 
This effect is illustrated in Fig. [5] (see the nodes in the 
gray boxes) and in Figs. [3j3 and c for i = 1. As can 
be visualised from Fig. [3J;, increasing the lattice position 
of the sink by one has the consequence that the connec- 
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tions between the new sink node and nodes (fc = 1, /) for 
1 < I < L ~ I (the first row in the rectangular subgraph 
of Fig. no longer belong to any s-t shortest paths. 
The argument can be easily generalized to any i. 




(2,2) 




1=0 1=1 1=2 



1=L-1 



1=L 



(b) 



15000^ 



5000" 




FIG. 6. (Colour online), (a) The number n(s, W(2,2) , i) of s-t 
shortest paths that cross node f (2,2) is determined from the 
product rule as the product of the number of s-t shortest paths 
in the two rectangular subgraphs highlighted in gray, (b) The 
number n(s, t) of s-t shortest paths between the source 

s and the sink t through the node {k, I) (where L — 10 and 
K = 12) has a saddle shape, with maxima at (fc = 1,Z = 1) 
and {k = K/2,l = L - 1). 



The number of s-t shortest paths passing through a set 
of given nodes is given by the product rule. The number 
of shortest paths connecting nodes s and t and passing 
through the node V(k,i) in the rectangular subgraph is 
given by: 



n{s, , t) = N{k, l)N{K/2 -k + l,L-l) 
k \\ f fK/2 -k + l 



I - 1 



L-l-l 



(11) 



where 1 < fc < K/2 and 1 < / < L - 1. Figure [6}i 
shows schematically how Eq. (fTTj) is derived, and Fig. [Sb 
is a plot of n(s, t) that illustrates that the con- 
centration of shortest paths is higher along the diago- 



nal of the rectangular layout. The number of short- 
est paths connecting nodes s = (fco,Zo) = (IjO) and 
t = (fcn+i,^n+i) = {K/2,L) and passing through nodes 
V{kuh),---,V(k^^i^) is given by: 

n 

= WN{ki+i-k, + l,li+i-k) 



i=0 



n 

i=0 



h+1 - fci + 1 
li+l — — 1 



(12) 



where 1 < fci < • • ■ < fc„ < K/2 and 1 < h < I2 < 
■ ■ ■ < In < L — \. From Eq. (fTTj) . the number of shortest 
paths containing edge 65 (^1) is defined by n{es,v(k i)) ~ 
"■(s, i'(fc,i), i) and is given by 



n{es,,^^J=N{K/2-k + l,L-l) 
'K/2 - fc + 1' 
L-2 



(13) 



Similarly, from Eq. (jl2p the number of s-t shortest paths 
passing through a set of nodes, we derive the num- 
ber of s-t shortest paths containing edges (^.^^ 1) and 
e(fe2,L-i),t that is defined as n(es,i,(fc^_i) , e^,^^^ t) = 
"(s,W(fe,i),W(fe2,L-i):*)> to be 



_i),t) =iV(fc2-fci + l,i-2) 

'k2 - fci + 1' 
L - 3 



(14) 



Equations pT|) to show that the number of s-t short- 
est paths passing through a node or a set of nodes can be 
written as a product of binomials. Each of these binomi- 
als is the number of s-t shortest paths when the s-t pair 
is located at the distance d™'" on networks with specific 
K and L values and is given by Eq. (fTO|) . 

To summarise, the number of s-t shortest paths de- 
pends on d and is given by a batch of sequences (see 
Fig. El), each having K/2 elements. If the s-t pair is 
at the smallest distance (F^™{K,L) given by Eq. ^ for 
a given shortest path length L, then the number of s- 
t shortest paths is maximal for that L. Increasing the 
distance between the s-t pair to d = d™™{K, L) + i (for 
1 < i < K/2 — 1) creates the same structure of s-t short- 
est paths as placing the s-t pair on a {K — 2i)-nearest 
neighbour network at the distance (i™™(if — 2i,L), which 
is the smallest distance at which the two nodes can be 
located on a {K — 2i)-nearest neighbour network so that 
the shortest path length between them is L. This im- 
plies that the problem of determining sink inflow on a 
isT-nearest neighbour network, as the distance between 
the s-t pair is varied, is reduced to the problem of calcu- 
lating sink inflow at specific distances d™"^{K,L) given 
by Eq. © as if and L are varied. Hence, we concentrate 
on determining the sink inflow rigorously for s-t pairs 
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placed at distances i)}|L=i,...,cxi on i^-nearest 

neighbour networks. From now on we will refer to a s-t 
pair at the shortest path length L to mean that the s-t 
pair is at the distance d = d''^''^{K, L), where L) 
is related to L in Eq. ([9]). As we will see in Sec. IIII CI 
the calculation of the sink inflow and of the path flows 
requires the number of s-t shortest paths in various sit- 
uations, as given by Eqs. (fTTj) to ^(T^ ■ 



C. Exact calculation of flows with the path 
counting methods in A'-nearest neighbour networks 
with one (s,t) pair 



contain shortest paths, 




]<j<i 



0<X<K/2-i+l-k, 



k=K/2-i+l 



o 

k=K/2-j+l 

o 

k=K/2 




1=0 



4 



1=1, 



1=1+1 



=L-1 1=L 



FIG. 7. (Colour online). Rectangular layout of the subgraph 
between the source (s) and the sink (t) nodes at iteration i 
of the MMF algorithm. The bottlenecks found at iterations 
1 ^ i < * a-re edges of s or t and are represented as thick 
black edges. The pair of bottlenecks found at iteration i of the 
MMF algorithm, are identified by the notation e^- ^+ (the 

western bottleneck) and e^- ^+ (the eastern bottleneck, not 
represented). When the bottlenecks are edges of the source 
or sink, oj^ = s and ef = t, respectively. 

Our goal is now to get analytical insights into the MMF 
allocation on A"-nearest neighbour networks with one s-t 
pair. In particular, we wish to determine all path flows 
as a function of the distance L between the source and 
sink nodes. The knowledge of all path flows also allows 
us to determine the incoming flow at the sink node as a 
function of L. The sink inflow can be determined as the 
sum of the path flows over the K/2 edges of the sink that 



F{K,L) = 



K/2 

1=1 pU,k)eP 



(15) 



where v^ ^ {K/2 - i + 1, L - 1), i = 1, . . .,K/2 are the 
K/2 neighbours of the sink that belong to s-t shortest 
paths. 

Such a problem becomes more interesting as we re- 
alise that the contribution to the sink inflow from each 
iteration of the MMF algorithm can be found by path 
counting methods alone. Indeed, the ratio of capacity to 
the number of paths for an edge on Eq. (IT|), is a function 
of both the edge capacity on the residual network and the 
number of unsaturated paths passing through the edge. 
In turn, the residual edge capacity depends exclusively 
on the number of paths that pass through the edge and 
are saturated at the previous iterations of the MMF al- 
gorithm. This implies that the calculation of sink inflow 
can be done rigorously from path counting alone, once 
we have found the pattern of the location of bottleneck 
edges at each iteration of the MMF algorithm. 

We start by assuming that all bottlenecks are found on 
edges of the source and sink nodes (see Fig. [T]). When 
K/2 is odd, this means that the row \K/4~\ in the rect- 
angular layout contains two bottleneck edges instead of 
one, and we observe a total of 2 [ii'/4] bottlenecks. The 
network is symmetric (i.e., the source and sink can be 
interchanged) and thus bottlenecks are positioned sym- 
metrically in pairs. We introduce the notation e^- ^+ 
and e^- ^+ , to refer to a symmetric pair of bottleneck 
edges found at iteration i. In general, we will refer to a 
bottleneck closer to the source e - + as a western bot- 

i ' i 

tleneck and a bottleneck closer to the sink e^- ^+ as an 
eastern bottleneck. Thus, at the i-th iteration of the 
MMF algorithm, the western bottleneck ^+ connects 
the source and node 

(16) 

and a symmetric eastern bottleneck e^- ^ will be located 
between the sink and node 

£7 ^ Vi^K/2^i+l,L~l}- (17) 

Since the bottleneck edges found at iteration i arc lo- 
cated on e^^+ and the symmetric edge e^- j, the path 

flow increment ([2]) at iteration i of the MMF algorithm 
is given by 



A/(') = min(/iW(e) 

eG-E 



)/ P^'^K, 



-) 



(18) 



for all c(*'(e^^+) 7^ 0. While the denominator of (fT8|) 
depends only on the number of unsaturated s-t short- 
est paths passing through ^+ at iteration i, the resid- 
ual capacity c*-*^ (e^, ^+ ) on the numerator is found recur- 
sively as a function of the path flow increment added to 
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unsaturated paths at iterations j < i of the MMF algo- 
rithm. This recursive calculation depends on the number 



V S.UjJ 



of s-t shortest paths passing through edge 
e + that have not been saturated until iteration j < i. 

1 i 

The path flow of each of these paths is increased by A/^-'^ 
at iteration j, and the capacity of all edges that these 
paths cross is reduced by Af^^\ 

Paths passing through ^+ that have been saturated 
in the previous iteration j < i of the MMF algorithm, 
have the property that each cross ^+ as well as one 
of the bottleneck edges e^- ^ found at a previous itera- 
tion q = 1, ... ,j — 1. For each q = 1, ... ,j — 1, there 
are n{e^ 



at iteration j crossing e 
Eqs. (Uni) and ^ as 



j) such paths given by Eq. p4|) . Hence, 

, which is the number of unsaturated paths 
, can be determined from 



i-i 



^N{K/2^{i+j) + 2,L-l), (19) 

where 1 < i + j — 1 < K/2 — 1, and from the paths 
counting formulas (fTO|) . (fT3|) and p4|) . and the location 
([T6|) and (fT7|) of the western and eastern bottlenecks, 



and 



J = Ar(q-i + l,L-2). 



(20) 



(21) 



To understand Eq. (|19p . observe that the s-i shortest 
paths crossing both ^+ and edges in row K/2 — j + 1 

of the rectangular subgraph (see Fig. [3]) also cross bot- 
tleneck e - f, and these paths are saturated at iteration 

j < iof the MMF algorithm. By symmetry, the s-t short- 
est paths crossing e^- ^ and edges in row j also cross the 
bottleneck e +, i.e. they are also saturated at iteration 

j. The effect of this pattern of saturated s-t shortest 
paths is that, for each iteration j < i, the MMF algo- 
rithm saturates two bottleneck edges, ^+ and e^- j, as 

well as all the paths that cross rows j and K/2 — j + 1 
of the rectangular subgraph. Thus, the effect of j itera- 
tions of the MMF algorithm is to remove 2j rows {j rows 
from the top and j rows from the bottom) of the original 
K/2-hy-L rectangular subgraph, creating an equivalent 
{K/2 — 2j)-hy-L rectangular subgraph. Note that the 
equivalent rectangular subgraph is only useful for path 
counting purposes, since the capacity of its edges will 
have been reduced after j steps of the MMF algorithm. 



Each of the 



unsaturated paths has its 



path flow incremented by A/^^'' at iteration j < i, ac- 
cording to Eq. So, each path in P'^^\e^ ^+) reduces 



the available capacity c*^^^(e^^+) by A/^-'^ according to 

Eq. ([3]). Since all edges have initial capacity c (see Eq. 
([7])), the path flow increment at iteration i of the MMF 
algorithm, given by Eq. (jlSp . can be written as a function 
of Eq. ([HI): 



A/ 














p^-^K^t) 





)A/(<?) 



Ello NiK/2 -{i + q)+2,L- l)A/(«) 



iV(A72-2(i- l),i- 1) 



(22) 



where A/^"^ ~ 0. Equation \22\ says that the incre- 
ment to path flows A/(*) at iteration i on a K/2-hY-L 
rectangular subgraph can be written as a recursive func- 
tion of the number of s-t shortest paths on a sequence of 
nearest-neighbour induced subgraphs with smaller node 
degree than the original rectangular subgraph. 
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FIG. 8. (Colour online). Plot of the ratio r^'-' between the 
path flow increment at iteration i and the path flow of paths 
that are saturated at iteration i, as a function of L for K = 24 
as the iterations i increase. The path flow is dominated by 
the path flow increment at the last iteration when r ~ 1. 

To gain insight into the Max-Min Fair path flows J*^* ' 
when all bottlenecks are edges of s or t, we plot in Fig. [5] 
the ratio 



A/ 



(0 



A/(') 



(23) 



for all paths that are saturated at iteration i. These 
paths that are saturated at iteration i are all the paths 
that pass through a bottleneck found at iteration i, or 
equivalently, all paths in P*^*^(eB). Thus, the ratio r^*^ 
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represents the relative contribution of the path flow in- 
crement at iteration i to the path flow assigned to paths 
that are saturated at iteration i. The observed r'*) — > f 
tells us that the path flows allocated by the MMF al- 
gorithm are dominated by the flow increment A/'*'' at 
the last iteration. A histogram of the path flows is 
thus well approximated by a plot of the number rj^"^^ = 
N{K/2 - 2{i - 1), X) - N{K/2 - 2i, L) of paths saturated 
at iteration i versus A/(*) oc l/N{K/2 - 2(i - f ), i - f ), 
and for L large A/*^*) ^ In other words, when 

the bottlenecks are edges of s or t and L is large, the 
histogram of path flows is a discrete set of points on a 
line with slope —1. 

Combining the number of s-t shortest paths and the 
path flow increment according to Eqs. (fTO|) and ([22]) . the 
MMF flow at the sink is then 



rif/41 



F{K, L)= J2 NiK/2 - 2{j - 1), L)A/ 



(i) 



(24) 



Combining the number of s-t shortest paths Eqs. ([T0|) and 
the number pO|) of paths passing through the western 
bottleneck, Eq. ([M]) can be simplified for large L to 



F{K,L) 



E 



K/2-2j+L\ 
L-1 / 



'K/2~2j+L-l\ 
L-2 ) 



(25) 



Using the asymptotic expansion of the binomial coeffi- 
cient, Eq. ([25]) can be simplified further to 



l)K/ 



2-2J + 1 



and thus 



lim F{K,L) = \K/A]c. 

L— >oo 



(26) 



(27) 



This asymptotic behaviour can be easily observed in Ta- 
ble |l] for small K. 



TABLE I. Sink inflow F{K, L) as a function of s-t shortest 
path length L for node degree 4 < K < 14, when all bottle- 
necks are edges of s or t. The asymptotic behaviour of the 
sink inflow is limL->oo F{K, L) = [Ay4]c. 

K F(K,L) 
4 
6 



10 
12 
14 



Lc 
L-1 

2c + 
2c + 
3c 
3c 

4c-|- 



(L-l)L'- 
L+4 „ 



2(l^ + 11L-2) 

L^+9L^+51L+H 
L^+Sf+SL^ -51,-6 



2(1, 
73 



-189L-'+88Z,-12) 



i:,<>+9L=+25L* + 15L-^-26L-^-24L 



So far, we have made the assumption that, at each 
iteration of the MMF algorithm, the algorithm saturates 



one edge of the source and another of the sink, up to 
a total of 2[A'/4] bottleneck edges. This assumption 
allows us to compute exactly the path flow increments 
(see Eq. ([H])) and the sink infiow (see Eq. (HH)). Knowing 
the pattern of the location of bottleneck edges gives us 
a clear advantage over solving the MMF allocation with 
the generic implementation of the algorithm described 
in Sec. IIIBI Indeed, the generic implementation of the 
MMF algorithm has to keep track of all paths and path 
fiows at each iteration. This process consumes increasing 
resources as K and L increase, owing to the number of s-t 
shortest paths growing polynomially or faster with either 
K or L. Thus, knowledge of the location of the bottleneck 
edges is particularly important for some regions of the 
parameter space, where the number of s-t shortest paths 
is large, and it is impractical to implement the MMF 
algorithm as described in Sec. Ill Bl 

Since our assumption has been that the bottlenecks 
are edges of s or t, it is then natural to investigate the 
regions of the parameter space (AT, L) where this assump- 
tion holds. To answer this question for each (A', L) pair, 
we compute the smallest iteration of the MMF algorithm 
such that the bottleneck edges are not edges of s or t. 
Since the bottleneck edges found at iteration i arc lo- 
cated on rows i and Ar/2 — i + 1, to find these bottlenecks 
we need to search for the minimum of Eq. ([1]) over all 
edges that link to a node in row i. This search can be 
greatly simplified when all 2{i — 1) bottlenecks found up 
to iteration i — \ are edges of s or t. Let us assume that 
the western bottleneck e - + , which is not an edge of s 
or t, connects the two nodes 



(28) 



where i < k2 < \K/A] , < A < A72-i + l-fc2 is an inte- 
ger, and the eastern bottleneck is defined in a similar way. 
When A = 0, we say that the bottleneck is horizontal be- 
cause it links two nodes that share the same row and are 
placed on adjacent columns on the rectangular subgraph 
in Fig. [3l When A 7^ 0, we say that the bottleneck is 
diagonal, because it connects two nodes that belong to 
adjacent columns, but different rows of the rectangular 
subgraph (see Fig. [7]). The search for bottleneck edges is 
simplified because we only need to search over horizontal 
edges. The proof consists in showing that the numerator 
(denominator) of Eq. ([1]) has a minimum (maximum) for 
horizontal edges, i.e. when A = and thus the search 
can be restricted to horizontal edges. 

We start by observing that the numerator of Eq. ([1]) 
is the edge capacity on the residual network at iteration 
i. To show that this edge capacity is the smallest when 
e^- ^+ is a horizontal edge, we demonstrate that it de- 
creases at iteration i by the largest value when e - + is 
horizontal, according to Eq. (|4]). Such decrease is given, 
at each iteration j < i oi the MMF algorithm, by mul- 
tiplying the path flow increment Eq. ^ by the number 
of paths that pass through e^- ^+ and are saturated at 
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that iteration j. The path flow increment is constant for 
all paths at each iteration. So we just need to show that, 
for all iterations j < i, the number of s-t shortest paths 
saturated at each iteration j, passing through e^- ^+ is 
larger when e - + is a horizontal, than when it is a di- 
agonal edge. Recall that here we are assuming that all 
bottlenecks found until iteration i — 1 are edges of s or t. 

Let us start with a K/2-hy-{L — 1) network. At itera- 
tion j < i, we have an equivalent {K/2 — 2{j — l)))-by- 
(L — 1) network. The paths that are saturated at iteration 
7 and cross e - + also cross i) the i-th western bot- 
tleneck es^„y and one of the edges ev^^^^_^^^ l-i),* f'^'" 
j < k < i; ii) the j-th eastern bottleneck ev^K/2^j+l L^i),t 
and one of the edges es,^,^ for j < k < i. Therefore, 
the total number of paths saturated at iteration j cross- 
ing e^-u,+ is 



at iteration j is found from Eqs. p2|) and ([281) as 



+ ^ "'(^S>t'(fc,l)lSw-,W+'^f(iC/2-j + l,I,-l) 

k=j+l 

K/2- 3 

"-(^s.^cj-.i) ' ^''(fc,I,-l).t)■ 



(29) 



The symmetry of paths from the source to the sink and 
vice- versa, together with Pascal's rule allow us to simplify 
Eq. dm to 



-)) 



j + f + ^ - 12 

L- 1-2-2, 
-J -l2 + k2-2\ /-j -l2 + k2-2 
h-2 ) \ h-l 

-j -l2 + k2-2 
I2-2 



L-h 



A 



-J + f + i - ^2 - A - /S2 - 2 

L-h -2 



-{L-h -3) 



(30) 



Observe that K/2 — j — k2 + 2 > 0. and thus the to- 
tal number of paths saturated at iteration j < i that 
cross e^- ^+ , which is given by Eq. (j30p . is maximum for 
A = 0, that is when e - + is a horizontal edge. This im- 
plies that the numerator of Eq. ([T]) has a minimum when 
the edge is horizontal. Next, we show that the denom- 
inator of Eq. ([1]) has a maximum for horizontal edges. 
To do this, we need to demonstrate that the number of 
unsaturated paths crossing e^- ^+ at iteration j < i is 
larger for horizontal edges than for diagonal ones. Note 
that all paths passing through the {K/2 — 2{j — l))-by- 
{L — 1) residual network are unsaturated at iteration j, so 
the number of unsaturated s-t shortest paths that cross 



unsat V'-'lj ■ .lj. 



(ea;r.c.+ ) 



k2 -j + l 

h-l 



K/2- j -k2- X + 2 
L-h~2 



(31) 



Equation (pT|) yields the maximum value for A = 0, 
i.e. for horizontal edges. Thus, the minimum of Eq. ([1]) 
is on a horizontal edge. This means that when the first 
2{i — 1) bottlenecks are edges of the source or sink, the 
2i-th bottleneck is either again an edge of s or t, or it is 
a horizontal bottleneck. 

Taken together, these results yield a simple procedure 
to locate the bottlenecks edges. Starting from the first 
row i = 1 on the rectangular subgraph, we search for the 
pair of bottleneck edges, i.e. the minima of Eq. ([T]) over 
the set of horizontal edges on that row. When i ~ 1, 
the bottlenecks are edges of s or t. Nevertheless, this 
is not necessarily the case when i > 1. In other words, 
there may exist an iteration i such that the bottlenecks 
previously found are edges of s or t, but the new pair of 
bottlenecks is not. We then know the exact location of 
the 2i bottleneck edges found up to iteration i. When i 
is the last iteration of the MMF algorithm, we can derive 
the Max-Min Fair fiows analytically from the path count- 
ing methods developed above, otherwise we still have an 
approximation of the MMF fiow that is exact up to it- 
eration i of the algorithm. Applying a similar reasoning 
to the one presented for the first inner bottleneck, it is 
possible to show analytically that all other bottlenecks, 
except the bottlenecks of the source (or sink), are hor- 
izontal edges as well when i) each row kj < ki has at 
least one bottleneck edge found in a previous iteration; 
ii) the last bottleneck to be found at each row kj < ki 
is an edge of the source (or the sink). The first condi- 
tion guarantees that, up to iteration i, all rows up to row 
ki have at least one bottleneck. The second condition 
guarantees that all s-t shortest paths crossing any edge 
in all rows kj < ki have been saturated, and thus the 
only unsaturated paths are in the {K/2 — 2{i — l))-by- 
(i — 1) rectangular subgraph. Since these two conditions 
simplify the path counting, we are able to compute the 
fiow distribution efficiently by a semi-analytical method. 
This method relies on an analytical analysis to restrict 
the search to horizontal links, combined with a numerical 
search over the same horizontal links. 



IV. RESULTS 

Taken together, the methods developed in Sec. Illll can 
now direct our investigations into Max-Min Fair (MMF) 
flows on nearest neighbour networks with one source and 
sink pair (s, t). We start by asking how sink inflow varies 
with node degree K, and the distance between the s and 
t nodes. As discussed in Sec. IIIIl we can measure the 
s-t distance in two different ways. On one hand, L is the 
shortest path length between s and t. On the other hand. 
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FIG. 9. (Colour online). Plot of the sink inflow as a function 
of (a) the s-t distance d and (b) the s-t shortest path length 
L for even K and 4 < 7^ < 14. The relation between d and 
L can be derived from Eqs. ([SJ and @. The solid symbols 
denote the sink inflow values common to both panels. The 
solid curves in (b) are computed assuming that all bottlenecks 
are edges of the source s or the sink t. The inset in panel (b) 
shows the approximate polynomial growth of the number of 
s-t shortest paths as L is varied for several K values. The sink 
inflow on a. s — t pair at the distance d = d'"*"(if = 14, L = 
2) +i for i = 1, ■ ■ ■ , 5 is the same as the sink inflow on a s — t 
pair at a smaller distance on a less dense network. The two 
cases are identified by circles of the same colour on panels (a) 
and (b). 



FIG. 10. (Colour online). Complementary cumulative distri- 
bution of the number of s-t shortest paths passing through 
each edge for L = 5, as d is varied. We considered only 
edges that are crossed by at least one path. The diversity of 
the number of paths crossing each edge is illustrated by the 
different distributions. When d = d"*"" (is', L) , the distribu- 
tion is broad and some edges are crossed by a large number 
of paths. The pattern of intersections among these paths 
constrains the solution of the MMF flow, because the paths 
share the capacity of the edges they cross. However, when 
d = d™'"(7i', L) -I- K/2 - 1, there is only one s-t shortest path. 
In this case, the MMF algorithm allocates the edge capacity 
to the path flow, because that path does not interact with 
any other. 



d{K, L) is the difference between node index labels for t 
and s, as illustrated in Fig. [51 when these labels are in- 
creasing from the source to the sink. The two distances 
are closely related, because there are K/2 s-t pairs at 
distances d = d'^'^iK, L), • • ■ , d'"'"(iC, L) -h K/2 - 1, but 
all of these s-t pairs are at the shortest path length L 
from the source s. Our first main result from Sec. IIII Bl 
is that we do not need to investigate the MMF algorithm 
for all s-t distances. In fact, when the s-t pair is at a 
distance d on a i^T-nearest neighbour network, the set of 
s-t shortest paths remains unchanged if the s-t pair is 
placed at a given smaller distance on a sparser network. 
To be more specific, displacing the s-t pair from the dis- 
tance d"'"'{K, L) to the distance d = d"""(/f, L)+i, for 
i = 1, • ■ • , K/2 — 1, is equivalent to placing the two nodes 
at the new distance d^™{K' ,L), on a A''-nearest neigh- 
bour network where K' /2 = K/2 — i. In other words, the 
s-t shortest paths are equivalent for the values of d on 
the {K, L) network and the values of K' on the sparser 
{K' , L) network, such that d-\- K' /2 is a constant for a 
given L. This result implies that the MMF path flows on 



if-nearest neighbour networks are uniquely determined 
by their values at the distances d^™{K^L), as K and 
L are varied. This effect can be observed in Fig. [9^, 
where we plot the sink inflow as a function of the dis- 
tance d between the source s and the sink t for several 
values of K. Such transformation from d™™{K,L) -\- i 
to c?"""(A'',i) explains the oscillatory pattern in the 
data points. To illustrate this oscillatory effect, we high- 
light with circles in Fig. [9h the sink inflow values for 
d = d"""(jr = 14, L = 2) -I- i, j = 1, • ■ • , 5, as well as 
the flow values for the corresponding d^^'°^{K',L = 2), 
demonstrating that the sink inflow is fully determined 
by its values at the distances L). 

Figure [9)3 is a plot of the sink inflow at distances 
d^™{K,L) as L is varied for several K values, where 
points common to Figs. [9^ and [Qb are represented by 
solid symbols. When all bottlenecks are edges of the 
source or the sink, the sink inflow can be computed with 
the procedure described in Sec. IIII CI and the result of 
such computation is given by the solid curves in Fig[3}3. 
Indeed, the computation of the path flow increments in 
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Eq. ([2]) , and thus the computation of the sink inflow, can 
be done by counting paths that pass through the bottle- 
neck edges, if the location of these bottlenecks is known, 
and this is our second main result (see Eqs. ([^^ and 
p4)) ). Furthermore, the path flow increments Eq. ([22]) 
can be determined in closed-form for each K on the re- 
gion of the {K, L) parameter space where the bottlenecks 
arc edges of s or t. 

Our choice of shortest path routing allows us to tune 
the level of interaction among the paths by varying the s- 
t distance d from d = o?™'" (maximum number of shortest 
path intersections) to d = d""" — K/2 + 1 (no intersec- 
tions). This can be observed in Figs. [2] and [51 which 
illustrate how the number of s-t shortest paths depends 
on d. It can also be observed in Fig. [TOl which shows the 
variation with d of the distribution of the number of s-t 
shortest paths passing through each edge. 

Whereas the path flows can be computed in close form 
if all bottlenecks are edges of s or t, it is still possible to 
find the path fiows from path counting methods alone, 
even when some bottlenecks are in less trivial locations. 
Indeed, our third main result is a semi-analytical proce- 
dure that reduces the computational complexity of the 
search for bottleneck edges. Brute force approaches to 
solving numerically the MMF flow assignment for large 
networks are beyond the capacity of the current genera- 
tion of computers. Indeed, for L large, the input N{K, L) 
to the MMF algorithm grows exponentially in size with 
K (for fixed L) and polynomially with L (for fixed K), 
as shown in the inset of Fig. [9)3. The search for bottle- 
neck edges on a generic network is done, at each iteration 
of the MMF algorithm, over all links of all nodes in the 
residual network, and that is a procedure of the order 
0(|1/| K). Our result simplifies the search procedure for 
nearest neighbour networks, so that its complexity be- 
comes 0(|F|), and therefore is independent of node de- 
gree K . The simplified search procedure consists in find- 
ing the minimum in Eq. ([2|) over the edges of the residual 
network that connect two nodes at distance d = K/2 (we 
say that these edges are horizontal). Once the location of 
horizontal bottleneck edges is known, the path counting 
methods developed in Sec. IIIII are the key to finding the 
path flows and the sink inflow. We find the location of 
the bottleneck edges, as well as the saturated paths, the 
path flows and the sink inflow, by repeating this proce- 
dure for each iteration of the MMF algorithm. Neverthe- 
less, the bottlenecks found are not necessarily edges of s 
or t, and we refer to such edges that do not link the source 
or sink nodes as inner bottleneck edges. The sink inflow 
data points in Fig. [SJd were computed with this procedure 
(also discussed in Sec. IIII C|) . and are thus exact results. 
The inset in Fig. [9)3 helps us to appreciate the power of 
the methodology: our method deals well with over 10® 
s-t shortest paths for K = \A and L = 40 easily, whereas 
a brute force version of the MMF algorithm as discussed 
in Sec. IIII would make heavy use of computational re- 
sources because all paths and path fiows must be stored 
during the execution of the algorithm. 




FIG. 11. (Colour online). Each cell corresponds to a {K,L) 
point in the parameter space, where K is even, and the source 
and sink pair are at a distance d™"{K, L) as given by Eq. ((9|. 
We partition the parameter space {K,L) in areas A^2i), such 
that after the first i iterations the MMF algorithm finds 2i 
bottlenecks that are edges of the source s or sink t for all cells 
inside an ^(2i) area. We use alternating red and blue coloured 
cells to distinguish neighbouring A(^2i) areas. The parameter 
space is partitioned into four regions. In region (a), delim- 
ited by the solid line, all bottlenecks are edges of the source 
s or sink t. Region (B), between the solid and dashed lines, 
is characterised by chains of inner bottlenecks followed by a 
bottleneck at s (or t) which still allows us to derive exact re- 
sults. Region ©, in lighter colours, is defined by the presence 
of gaps in the row number of consecutive bottleneck edges. 
Region @ is a subset of © where the gap occurs between a 
bottleneck of s or t and the first inner bottleneck. 



The case for our procedure becomes even more com- 
pelling when we investigate the structure of the parame- 
ter space {K, L). To show this, we construct a parameter 
space diagram, where we associate with each (K, L) pair 
the number of bottlenecks that are edges of s or t up 
to the iteration when the first inner bottleneck is found. 
The result, which we plot in Fig. [Til is a partition of the 
parameter space in areas A2i, such that for K and L val- 
ues in there are 2i bottlenecks which are edges of s 
or t. For L sufficiently large, the behaviour of the MMF 
algorithm is described by the area A2. In other words, 
the bottlenecks are edges of s or t only for the first iter- 
ation of the MMF algorithm. This means that the two 
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bottlenecks found at the first iteration are edges of s or 
t, but the bottlenecks found at the second iteration are 
already inner bottlenecks. In fact, when L is large, we 
only understand well the first iteration of the MMF al- 
gorithm, and this opens the possibility that the pattern 
of the location of bottlenecks in the area A2 of the pa- 
rameter space may be much more complicated than we 
have been able to describe. Four regions of the parameter 
space illustrate what we know about the MMF algorithm 
in nearest neighbour networks. The first region (a), de- 
limited in Fig. [H] by a black line, is the set of points in 
the parameter space {K, L) such that all bottlenecks are 
edges of s or and this is where we are able to derive 
analytical results. The second region @, located between 
the solid and dashed lines, limits the {K, L) values where 
the semi-analytical methods developed in Sec. IIII CI can 
solve all iterations of the MMF algorithm. These semi- 
analytical methods enable the analysis of the path flows 
and the sink inflow as L grows. Our main flnding is that 
the calculation of sink inflow for L large in region @ (see 
Eq. (|27p ) is an upper bound to the sink inflow in region 
(b), as can be observed in Fig. [5)3. This can be explained 
by the pattern of the location of bottleneck edges in (B) 
(see Fig. [TT|) . More precisely, the capacity constraints 
are larger in (b) than in (a), because the pattern of bot- 
tlenecks in (b) contains the pattern of bottlenecks in (a). 
The third region © is characterised by the presence of a 
gap in the row number of two consecutive bottlenecks on 
the rectangular subgraph. Region @ is a subset of region 
© where the row gap occurs between a bottleneck adja- 
cent to s or t. This 'gap' in the row number of bottleneck 
edges implies that a subsequent bottleneck has to appear 
in the rows that have been skipped. Such detailed knowl- 
edge of the location of bottleneck edges in the parameter 
space of nearest neighbour networks is one of the main 
contributions of this paper, but it also highlights the lim- 
its to our current knowledge of the MMF algorithm, as 
our methodology is not able to continue computing path 
flows efficiently for iterations where row gaps occur. The 
main limitation is that the implementation of the path 
counting method becomes extensively complicated with 
the appearance of row gaps. 

A natural question to ask is then: can network topol- 
ogy influence the Max-Min Fair allocation of path flows? 
To gain insight into this problem, we plot the histograms 
of path flows in Fig. [T^ for several {K, L) points inside 
the regions of the parameter space where we can solve the 
MMF algorithm exactly (regions © and © in Fig. [12]). In 
the region @ of parameter space where all bottlenecks are 
edges of s or the histograms show that the frequency 
of path flows decays approximately as a power-law with 
exponent —1 (sec Fig. 112b ,). Two mechanisms contribute 
to the way that fair flows are allocated in this area. First, 
the path flows are dominated by the path flow increment 
at the last iteration, when the paths are saturated. From 
Eq. we can write 

fi'*} ^ ^fi^n ^ i/N{K/2 - 2{i* - 1),L). (32) 
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FIG. 12. (Colour online). Histograms of the path flow after 
the MMF algorithm terminates, (a) When the bottlenecks 
are edges of either the source or sink nodes (or we observe 
only a small number of inner bottlenecks) the histogram de- 
cays approximately as a power law. (b) As the number of 
inner bottlenecks grows, we observe a scattered distribution 
that is caused by different path flows, specific to each in- 
ner bottleneck. Two bottlenecks that are on the same row 
of the rectangular subgraph have similar path flows. As a 
consequence, path flows of bottlenecks on the same row are 
clustered on the histogram. 



Second, the number r](i) of s-t shortest paths that are 
saturated at each iteration is of the order of magnitude 
of the number of all paths in the residual network, that 
is, 
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~iV(X/2-2(r -1),L). 



(33) 



Hence, at each iteration i*, the path flows, given by 
Eq. (|32|) . are approximately inversely proportional to the 
number of saturated paths at iteration i*, determined 
by Eq. In other words, ~ l/v^''\ and the 

histogram of path flows decays as a power-law with ex- 
ponent — 1, as seen in Fig. [T2k. On the other hand, the 
presence of inner bottlenecks in the region ® of the pa- 
rameter space introduces a deviation from the power-law 
decay (see Fig. [T2b). Indeed, paths crossing the bottle- 
necks that share the same row of the rectangular sub- 
graph have similar path flows, and thus the frequencies 
of paths flows are clustered on the histogram. Taken 
together, these flndings show that power-law allocations 
can be fair, even when network capacity is uniform, which 
is a counter-intuitive and unexpected result. 
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CONCLUSIONS 



We analysed the Max-Min Fair flow allocation algo- 
rithm on nearest neighbour networks with uniform capac- 
ity and node degree K. Our study focused on one single 
source and sink pair, placed at shortest path distance L. 
Transport between the source and the sink can only oc- 
cur along the shortest paths, each of which is assigned 
a path flow by the Max-Min Fair algorithm. We note 
that we have analysed the {K, L) parameter space under 
the assumption that, for each L value, the s-t pair is at 
the distance d = d^''^{K,L) as given by Eq. ®. Other 
distances d = d"'"\K,L) -f- i (for 1 < i < K/2 - 1) can 
be mapped to an s-t pair located at a shorter distance 
on a sparser network (as discussed in Sec. IIII B|l . 

Such configuration led us to analyse how the tuning of 
the parameters, K and L, affects the network through- 
put between the source and sink pair, as well as the fair 
allocation of individual path flows. Figure [9] shows the 
behaviour of the network throughput as a function of the 
distance between the source and the sink. 

We found that the parameter space {K, L) can be par- 
titioned into areas that correspond to distinct patterns in 
the location of congested edges, i.e. bottlenecks. When 
the bottlenecks are edges of the source s or sink i, we de- 
rived a recursive relation that yields the sink inflow and 
the path flows analytically. Furthermore, in the limit 
L — ^ oo, the sink inflow grows with node degree as the 
piecewise constant function \\vaL^aoF{K,L) = \K/A\c. 
This result is exact for node degree K < A, when all 
bottlenecks are edges of the source or the sink, and we 
showed that [i<r/4]c is an upper bound to the sink inflow 
when K > A and L is large. In contrast, the maximum 
flow that can be allocated between the source and the 
sink node pair is Kc/2 [49| . Hence, in the limit i — >■ oo 
there is a decrease in throughput larger than 50% of max- 
flow, and this limitation of the Max-Min Fair algorithm 
is well-known in the literature (see e.g. [H, H^l). Indeed, 
Betsimas et al. [l^ established an upper bound for the 
loss of throughput in a very general setup, and when ap- 
plied to our case with L — > oo these results imply 100% 
reduction of the throughput. Moreover, numerical stud- 
ies show that when the number of s-t pairs is large and all 
paths are allowed, the result is a reduction from max-flow 
comparable to the case of fixed shortest paths [s^ . How- 
ever, it should be noted that both our results and those 
of Nace et al. [ssj are found when paths have many inter- 
sections and, hence, are not valid if the network operator 
would choose the paths not to intersect. Our results can 
thus be seen as going further, by providing an exact de- 
scription of the flow reduction in specific areas of the 
parameter space. Taken together, these findings suggest 
that implementations of the Max-Min Fair algorithm in 
very large real world networks should be restricted to 
situations where throughput can be sacrificed to achieve 
fair allocations. 



For small values of L, all the bottlenecks are edges 
of s or < and path flows can be computed analytically. 
As the distance L between the source and the sink in- 
creases for a fixed K ^ the pattern of intersections among 
shortest paths changes and the location of the bottle- 
necks becomes less regular. To address the new patterns 
of the location of the bottlenecks, we then derived a semi- 
analytical method that combines analytical results that 
reduce the complexity of the search for bottlenecks, with 
numerical computations for the search. As a result, we 
were able to solve the Max-Min Fair allocation for large 
nearest neighbour networks when K < despite the 
very large number of shortest paths involved in the com- 
putations. 

When all bottlenecks arc edges of the source or the 
sink, we find that the histogram of path flows decays 
approximately as a power- law with exponent —1. This 
uneven distribution is unforeseen, because the Max-Min 
Fair algorithm assigns path flows fairly, i. e. the least well 
off get as much as possible under the capacity constraints. 
This counter-intuitive result that power-law allocations 
can be fair is a consequence of the constraints placed on 
the Max-Min Fair allocation by the combination of the 
network topology and the routing over shortest paths. 
Moreover, we observed deviations from this power-law 
decay in regions of the parameter space where some bot- 
tlenecks are not edges of the source or sink. These unex- 
pected results suggest that network designers should be 
aware of the interplay between the structure of transport 
routes and network topology, which is a crucial factor in 
the fair allocation of network flows. 

Finally, the analysis developed in this paper may be 
extended in several directions. A possible line of enquiry 
is to study the network throughput as a function of the 
number of distinct sink-source pairs. Still another possi- 
ble step is to investigate the effect of network topology 
on Max-Min Fair flows. Further studies arc required to 
understand whether the reduction in throughput that we 
have found is specific to the choice of fixed shortest paths, 
or whether it is a more general result. Our study is just a 
first step towards a rigorous understanding of fair alloca- 
tions in complex networks. Due to the importance of the 
topic, further insights could be reached by analysing ran- 
dom graphs or real world network topologies. It would 
also be important to gain a better understanding of the 
relation between different routing strategies and the net- 
work throughput. Last, but not least, it would be worth- 
while to consider other fair allocation schemes, such as 



the Proportional Fairness [24 
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